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Abstract 

The usual notion of separability has to be reconsidered when applied to states describ- 
ing identical particles. A definition of separability not related to any a priori Hilbert 
space tensor product structure is needed: this can be given in terms of commuting 
subalgebras of observables. Accordingly, the results concerning the use of the quantum 
Fisher information in quantum metrology are generalized and physically reinterpreted. 



1 Introduction 

Quantum interferometry is one of the most useful measurement technique used in quantum 
metrology: it exploits the high sensitivity of certain quantum systems to small changes of 
external conditions, due to their intrinsic quantum coherence. Although various implemen- 
tations of this general measurement procedure have been realized, the aim is always that of 
estimating the relative phase of the two "arms" of the interferometer and great efforts have 
been pursued to enhance the sensitivity of this phase determination [1, 2J. 

The accuracy of the whole estimation procedure [3|-[n] depends on the number of resources 
at disposal, i.e. on the total number N of particles involved in the interferometric measure. 
It turns out that the minimum uncertainty in the phase estimation that can be achieved by 
feeding the interferometer with "classical" states can not exceed the shot noise limit, also 
called the standard quantum limit, where the precision scales as l/y/N. 

Various strategies have been devised in order to beat this limit, typically adopting non- 
standard measurement protocols that make use of quantum correlated (entangled), non- 
classical statesQ In this way, at least in line of principle, the so called Heisenberg limit in 



1 The literature on the subject is vast, dealing both with optical as well as atom interferometers; for a 
partial list, see [7]- [33] and references therein. 
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accuracy can be reached, where the uncertainty in the phase estimation scales as 1/N rather 
than l/VN. 

A particularly promising implementation of an interferometer that could reach such a high 
sensitivity is that of a Bose-Einstein condensate (BEC) trapped in a double- well 

optical potential. Although similar to more traditional constructions that use beams of 
particles travelling along separate paths, this new realization of the interferometry paradigm 
allows a precise control on the preparation and dynamics of the system. 

The clear advantage in using BEC-based interferometers is that all particles in the appa- 
ratus, being in a condensed phase, share the same quantum state; this fact, together with 
the possibility of preparing them in a suitably entangled state, have been shown to lead to 
the possibility of phase estimation with sub-shot-noise accuracy [22]- [32]. However, these 
particles, being spatially confined bosons, behave as identical particles, a fact that seems to 
have not been properly taken into account in the derivation of the above mentioned results. 

Purpose of the present investigation is to study the effects of particle indistinguishability 
in relation to the metrological use of BEC-based interferometers for phase estimation. Since 
entanglement appears to be a fundamental resource in reaching sub-shot-noise accuracies, one 
first needs to generalize this notion to states describing identical particles. This is presented 
in Section 2; there, we shall see that one needs to abandon the usual definition of separability 
expressed in terms of tensor product of one-particle Hilbert spaces, in favor of a more general 
definition based on the choice of commuting algebras of observables. In Section 3 we shall 
explicitly discuss how this generalized notion of entanglement, based on a dual description in 
terms of operators rather than states, works in the case of a BEC confined in a double-well 
trap. Finally, in the light of this discussion, in Section 4 we shall critically re-examine the 
theoretical results concerning the use of the notion of quantum Fisher information in getting 
sub-shot-noise accuracies in quantum metrological phase estimation. The final Section 5 
contains further discussions and remarks. 



2 Identical particles entanglement 

The usual notion of entanglement for states of a system of N distinguishable particles makes 
use of the natural tensor product structure of the iV-body system. The total Hilbert space 
% of the system can in fact be decomposed as H = Hi <S> Hi ® • • • ® Hn in terms of the 
various single-particle Hilbert spaces Hi, i = 1, 2 . . . , N. Then a state for the iV-body system, 
represented by a density matrix p acting on H, is said to be separable if it can be written 
as a convex combination of single-particle states, 

p= ^Pkp k 1] ®p k 2) ® •••®pf ) , p k >o, ^2 Pk = i, (i) 
k k 

where p^ represents a density matrix for the z'-th particle; otherwise, the state p for the N 
distinguishable particles is said to be entangled. 
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The tensor product structure used in ([I]) reflects on one hand the natural multi-partition 
of the system into its N parts, and on the other the possibility of identifying linear local 
operations on the system, i. e. maps of the fornjl 

A = A« <g> A® ® . . . ® A<^ , (2) 

where the linear map : i— > A^ jp^l affects only the i-th particle density matrix. 

Coming now to the case of a system of N identical particles, one easily realizes that states 
of the form ([T]) are not allowed quantum states for the system. Indeed, if one could assign a 
density matrix pM to the first particle, p^ to the second one and so on, we would have a 
way to distinguish them by means of their identifying states. 

On the other hand, recalling the standard rules of quantum mechanics, a pure state \%p) of 
N identical particles must be a symmetric or antisymmetric combination of tensor products 
of A^-single particle vector states, while a mixed state, i.e. a density matrix, must be a 
linear convex combination of projections \ip)(ip\ onto such symmetrized or antisymmetrized 
vectors gDH2]ll 

In particular, this result excludes the possibility of curing the intrinsic distinguishability of 
the particles when described by states as in (JTJ), by symmetrizing or antisymmetrizing them: 
indeed, in the end, this would always be a mixture of states where particles can be effectively 
distinguished. In other words, density matrices describing identical particles must be made 
of projections onto pure states of the appropriate (symmetrized or antisymmetrized) Hilbert 
space. As a consequence, not even the so-called "symmetric states" of the form 

P = ^Pk Pk ® Pk ® • • • ® Pk , (3) 

k 

where, unlike in ([I]), all particles in each term of the convex combination are described by the 
same one-particle density matrix, are in general admissible states for a system of identical 
particles. 

The example of a system of two qubits can further clarify these points. The Hilbert space 
of two distinguishable qubits is C 2 <g> C 2 = C 4 , and therefore four-dimensional. Instead, the 
Hilbert space for two identical qubits is a symmetric three-dimensional subspace of C 4 in 
the case of bosons, or an antisymmetric one-dimensional subspace for fermions. Given an 
orthonormal basis in C 4 , |+,+), |+, — ), | — ,+), |— , — ), where + and — indicates the two 
single-qubit levels, the symmetric sector is spanned by |+,+), |+,+), and the combination 
(|+, — ) + | — , +))/v / 2, while the antisymmetric one is generated by (|+, — ) — | — , +))/y/2. It 
is now easy to realize that density matrices of the type ([I]), p = ^2 k PkPk^ ® P k 2 \ can not 
describe states of two identical qubits. Indeed, a qubit state pf 1 is represented by a 2 x 2 

2 In order to be physically meaningful, these linear transformations must be realized by trace-preserving, 
completely positive maps |38l 139] ; only in this case, normalization, Tr[p] = 1, and positivity, p > 0, of the 
states p of the TV-body system are assured in any situation (for a discussion, see [10]). 

3 Similarly, also the observables must be symmetric or antisymmetric under the permutation of all particles; 
in particular, collective, TV-body operators are surely admissible observables (see also the note on page 9). 
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density matrix, so that any element p k (g> p k in the previous convex sum is a 4 x 4 matrix, 
as is the resulting combination p. As such, it can not be written in general as a convex 
combination of solely projections onto symmetric states or in terms of the projection onto 
the antisymmetric one, as it would be instead required for density matrices describing two 
identical qubits. Clearly, the same conclusion applies also for two qubit states of the form 

Y,kPk Pk®Pk- 

Given that the notion of separability based on the decomposition ([1]) is no longer applica- 
ble in the case of a system of identical particles, one is forced to generalize it. This problem 
has already been addressed and discussed in the literature, although in different physical 
contexts, e.g. see [I3]-[56j. The way out seems to reside in a "dual" point of view, where 
emphasis is given to the algebra of observables, instead of the set of states; the connection 
between these two "points of view" is given by the expectation value map that allows to 
express the average of an observable O as its trace with the corresponding density matrix, 
(O) = Tr[pO]. We shall first give an abstract definition of this generalized notion of separa- 
bility and then, in the next Section, apply it to the specific case of condensed bosons trapped 
in a double-well potential. 

Let us consider a system of particles whose pure states span the Hilbert space H and 
denote by B{%) the algebra of all bounded operators on it; the observables of the system 
clearly belong to this algebra. We shall introduce the notion of (bipartite) separability by 
considering couple of commuting subalgebras of B(7i) instead of focusing on partitions of 
the system quantum states. We then introduce the following definitions: 

• An algebraic bipartition of the algebra B(H) is any pair {A\ : A-i) of commuting subal- 
gebras of B{H)n 

• An element (operator) of B{H) is said to be local with respect to the bipartition 
(.Ai, .A2) if it is the product A\A 2 of an element A x of A\ and another A 2 of A 2 . 

• A state u on the algebra B(H) will be called separable with respect to the bipartition 
(Ai,A2) if the expectation u(A 1 A 2 ) of any local operator A^A 2 can be decomposed 
into a linear convex combination of products of expectations: 

co(A 1 A 2 ) = J2^co i k 1) (A 1 )J k 2 \A 2 ) , A fc >0, ]TA fc = l, (4) 

k k 

where u k and u> k are states on £>("H); otherwise the state u is said to be entangled 
with respect the bipartition (Ai,A 2 ). 

As we shall see in the following, this generalized definition of separability results meaningful 
both for systems composed of distinguishable particles and for those made of identical ones. 

4 Notice that the two subalgebras Ai and A2 need not reproduce the whole algebra B(H), i.e. in general 
Ai U Ai <zB(H). In this respect, the term "bipartition" is not strictly appropriate and has been adopted 
for sake of simplicity. However, in the case of the system discussed below, the considered mode partitions 
actually generate the whole algebra B(H). 
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In particular, in the case of two finite-level systems, as the two qubits example discussed 
before, the above definition reduces to the usual notion of separability. Indeed, it is enough 
to choose the subalgebras A\ and A2 to coincide with the 2x2 matrix algebras of the single- 
qubits and take as operation of expectation the usual trace operator over the corresponding 
density matrix, u(AiA 2 ) = Tr[pAi <g> A 2 ]] clearly, this expectation can be written as in (TJJ 
only if the system state p is as in ([1]) with N = 2. Finally, generalization of the above 
definition to the case of more than two partitions is straightforward. 



3 Identical bosons in a double-well trap 

We shall henceforth focus on the study of a system of ultracold bosonic atoms trapped in an 
optical lattice. As mentioned in the Introduction, the high degree of accuracy and sophistica- 
tion reached in manipulating such systems allow very accurate interferometric measurements, 
making them a very attractive set-up for exploring and testing new paradigms in quantum 
metrology and many-body physics [57] 158]. 

The dynamics of cold atoms confined in an optical trap can be very well approximated by 
a Bose-Hubbard type Hamiltonian [SHIED]; in the case of a double- well potential, it reads: 



Hbh — + £2 a \ a 2 + U (a\ai) 2 + (c4a 2 ) 2 — J (a\a 2 + aia\) , (5) 

where 0^2, a\ 2 annihilate and create atom states in the first, second well, respectively, and 
satisfy the Bose commutation relations [a*, aj] = 5ij, i,j = 1,2. Of the contributions to 
Hbh, the last one corresponds to a hopping term depending on the tunneling amplitude J, 
the first two are due to the trapping potential and are proportional to the depth e 12 of the 
wells. Finally, the third term, quadratic in the number operators ajaj, takes into account 
repulsive Coulomb interactions inside each well. 

We are using here a formalism of second-quantization, where the creation operator fo(x) 
of an atom at position x can be decomposed in general as 

00 

i=l 

and a\ creates an atom in the state \<pi) = a]|0) with wavefunction <pi(x) = {x\(j>i)\ the set 
{l^i)}"^ is therefore a complete set of orthonormal single-particle atom states. The Bose- 
Hubbard Hamiltonian flS]) results from a tight binding approximation, where only the first 
two of the basis vector are relevant; in this case (fii^ix) are orthogonal functions, <pi localized 
within the first well, 02 within the second one. 

The total number N of atoms is conserved by ([5]). Therefore, the Hilbert space of the 
system is iV + 1-dimensional and can be spanned by Fock states of the form 

\k,N-k):=-%lj&=\0), 0<k<N, (7) 
V k\^/(N — k)\ 
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describing the situation in which the first well is filled with k atoms, while the other one 
contains N — k particles. They are obtained by the action of the creation operators on the 
vacuum. These states turn out to be eigenstates of the Bose-Hubbard Hamiltonian when 
the tunneling term can be neglected. 

Notice that in this formalism, symmetrization of the elements of the Hilbert space, as 
required by the identity of the particles filling the two wells, is automatically guaranteed by 
the commutativity of the two creation operators |41j . Furthermore, all polynomials in ax, a[ 
and similarly all polynomials in a 2 , a\ (together with their respective norm closures), form 
two algebras, Ax, respectively A 2 , that result commuting: [Ax, A2] = 0; they are subalgebras 
of the algebra A of all operators on the Fock space spanned by the states ((7j)Jl According to 
the definition introduced at the end of the previous Section, they define a bipartition of A 
and therefore can be used to provide the notion of separability for the states describing the 
identical atoms in the trap, generalizing the one usually adopted for distinguishable particles. 

With respect to this natural mode bipartition, (Ax, A2), the Fock states (JTj) turn out to be 
separable states. Indeed, for any polynomial operator Ax(ax,a\) G Ax and A 2 (a 2 ,a\) G A 2 , 
the expectation value of the product A X A 2 on such states can always be written as the 
product of the separate averages of Ax and of A 2 ; explicitly, one has: 

(k,N-k\AxA 2 \k,N-k) = (0\a k A 1 (a\) k \0) (0\a^ k A 2 (4)^0) 

= (k\Ax\k) (N — k\A 2 \N — k) , (8) 

where \k) := (a\) k /Vk\\0) and \N—k) := (a\) N ~ k / \J(N — k)\\0) are single-mode Fock states. 

As a consequence, separable with respect to the bipartition (AijAz) are also those mixed 
states that are diagonal with respect to the Fock basis (J7|), i.e. density matrices of the form: 

TV N 

p = ^p k \k,N -k)(k,N -k\ , Pk >0, £}p fc = l. (9) 

k=0 k=0 

Actually, all states separable with respect to the bipartition (Ax,A2) must be in diagonal 
form with respect to the Fock basis. An abstract proof of this fact can be easily given. In 
fact, one first observes that the algebra generated by the polynomials in the operators ax, a[ 
and a 2 , a\ coincides with the whole algebra of operators A. Then, being pure projections, 
the states \k, N — k)(k, N — k\ are the only extremal ones in the convex set of states satisfying 
the separability condition (j3j); thus, all others are necessarily a convex combination of them. 

Nevertheless, a more direct proof can also be explicitly worked out. A generic density 
matrix representing a state for a two-mode system of N identical bosons is of the form: 

N N 

P = ^ Pki\k,N -k)(l,N -l\ , p kk >0, ^ Pkk = l. (10) 

k,l=0 k=0 



5 Strictly speaking, polynomials in bosonic creation and annihilation operators are not bounded, and the 
so-called Weyl algebras of the corresponding exponentials should be used; however, the following discussion 
is not affected by working with algebras of unbounded operators. 
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Assume now the state p to be separable with respect to the bipartition (Ai, A2); this means 
that for any local operator A1A2, with A\ G A\ and A2 G A2, one must have: 



Ti[pA 1 A 2 }=Y,^np^A 1 }Tr[p^A 2 } , 



(11) 



where p 



'm j 



1, 2, are density matrices of the generic form ([TO]) . Choose now Ax = (a[) m a", 



with m < n and = (c4) r a 2' with s < r and m + r = n + s; then, each factor in the sum 
above identically vanish. Indeed, one finds: 

JV TV 

Tr[p« A t ] = £ p£2w < fc > * - # - = J>2 fc * (fc^l*) = , (12) 



k,l=0 



k=0 



and similarly, Tr[p^A 2 ] = 0. However, using (Til)]) , one also have: 



/v 



iV-s 



Tr[p AiA 2 ] = P kl (^l^ilO " ^2|iV - /) = ^ Pfc, fc -n+ m a m> „(A;) /3 r , fl (fc) 

k,l=0 k=n 

where the coefficients a and /3 are explicitly given by: 

a m ,n( k ) = H k - 1) • ■ ■ (k ~ n + 1) ■ (k - n + l)(k - n + 2) ■ ■ • (k - n + m) 
(3 ryS (k) = \(N - k)(N - k - 1) ■ • • (JV - k - s + 1) 

x(iV - k - s + 1)(JV -k-s + 2)---(N -k-s + r) 



(13) 



1/2 



1/2 



Since Tr[p vanishes by hypothesis, by choosing n = N—s, one finally gets the condition: 



m\n\ (JV - m)\ (JV - n)\ 



1/2 







m < n 



(14) 



a similar condition holds for the case m > n. As a consequence, p mn = for all m + n, and 
thus the original density matrix in f llOp results diagonal in the Fock representation]^ 

However, most observables of physical interest are non-local with respect to the bipartition 
(-4.1, ^4.2) given above. As we shall see in the next Section, the approach to phase estima- 
tion with BEC-based double-well interferometers is mainly based on the following collective 
bilinear operators: 

1 1 1 

J* = 2 ( a i fl2 + aia 2) ' J v = 2~- ( a i a 2 _ a i4) > Jz = 2 ( a i a i ~~ a 2 a 2) > (I 5 ) 



satisfying the su{2) algebraic relations: 

\_Jxi Jy\ iJz 1 \pyi J z\ i"Jx 1 \_Jzi Jx\ iJ't 



(16) 



6 Furthermore, one easily checks that the states © satisfy the separability criterion for continuous variable 
systems introduced in (|61j). 
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The effect of phase accumulation inside the interferometer can in fact be modeled through 
the action of the exponential of a linear combination of these operators. It is interesting to 
notice that, although the operators in (1 15]). being linear combinations of both ai, a[ and a 2 , 
a\, are clearly non-local with respect to the bipartition (Ai,A2), their exponentials are not 
all so. While e %dJx and e iejy , 9 G [0,27r], are non-local, the exponential of J z turns out to be 
local: 

e «6>J z _ ^,8a\a\/2 _ ^-i0a\a 2 /2 ^i9a\a 1 /2 ^ j^^ ^-ida\a 2 /2 ^ (^-7) 

Further, by changing the bipartition of the algebra of all operators A, one can transfer 
this property to another of the operators in (I15p . For instance, let us introduce a new 
set of creation and annihilation operators b\, bi, i = 1,2 through the following Bogolubov 
transformation: 

a 1 + a 2 , ai - a 2 /1oX 

and their hermitian conjugates. Correspondingly, the three operators in ( TT5]) can be equiva- 
lently rewritten as: 

4 = \(b\bi - b\b 2 ) , J y = i(6 x «4 - &l& 2 ) , J 2 = + b\b 2 ) . (19) 

In analogy with what has been done before, one can now define a bipartition (Bi,B 2 ) of 
the full algebra A using the mode operators b\, bi instead of a\, ai. In this case, it is the 
exponential of the operator J x that now turns out to be local: 

e iej x = e ieb{ bl . e -ieb\b 2 ^ e ieb\b x e ^ ^ e -m\b 2 e ^ _ 

This explicitly shows that an operator, local with respect to a given bipartition, can result 
non-local if a different algebraic bipartition is chosen. 

From the point of view of the states, the above Bogolubov transformation corresponds 
to a change of basis in the Hilbert space, from the one consisting of spatially localized 
states (a\ and create and destroy particles in the two wells), to the one spanned by 
their (spatially non-local) superpositions, as, for instance, &{|0) = [a\\0) + (4|0)] /\/2 and 
&t| ) = [aJlO) - 4|0)]/V2, which are energy eigenstates of the Bose-Hubbard Hamiltonian 
in the limit of a highly penetrable barrier. As a consequence, the Fock states in ([7]) result 
entangled with respect to this new bipartition (Bi,B 2 ); indeed, one finds: 

|^-*> - ^-jmTwtt C" y- 1 ^-^ ®>"-~ |0> ■ (21) 

so that \k, N — k) is a combination of {B\, i32)-separable states. A similar conclusion applies 
to the mixed states ([9]). This explicitly shows that in the case of systems of identical particles, 
the notion of separability can not be given abstractly, but must be associated to a specific 
bipartition of the full algebra A of operators. Because of their physical meaning, in the 
following we shall refer to the (Bi, B 2 ) bipartition as the energy bipartition, and call spatial 
bipartition the original {A\, A 2 ) one. We shall now apply these considerations to the problem 
of phase estimation in quantum interferometry with systems of identical particles. 
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4 Quantum metrology with identical particles 



Ultracold atoms trapped in a double-well optical potential realize a very accurate interfer- 
ometric device: indeed, state preparation and beam splitting can be precisely achieved by 
tuning the interatomic interaction and by acting on the height of the potential barrier. The 
combination of standard Mach-Zhender type interferometric operations, i.e. state prepa- 
ration, beam splitting, phase shift and subsequent beam recombination, can be effectively 
described as a suitable "rotation" of the initial state p m by a unitary transformation [91 [TTj: 

Pin ^ Pe = U e p in Ul , U e = e w J " . (22) 

The phase change is induced by the generator 

J n = n x J x + n y J y + n z J z , n\ + n 2 y + n\ = 1 , (23) 

i.e. by a general combination of the collective operators introduced in ffl5|) . In practice, 
the state transformation p in y p g inside the interferometer can be effectively modeled as 
a pseudo-spin rotation along the unit vector n = (n x ,n y ,n z ), whose choice depends on 
the specific realization of the interferometric apparatus and of the adopted measurement 
procedure. 

General quantum estimation theory allow a precise determination of the accuracy on which 
the phase change can be measured in such devices: the accuracy AO with which the phase 
6 can be obtained in a measurement involving the operator J n and the initial state p m is 
limited by the following inequality [3] -[6]: 

A9 > - 1 , (24) 
- >F[p in , J n ] 



where the quantity F[p- m , J n ] is the so-called quantum Fisher information. One can show 
that in general [62J 

F[p in , J n ] <4(A Pin J n ) 2 , (25) 

1 /2 

where A Pin J n = [(Jn) — (J n ) 2 ] is the variance of the operator J n in the state pi n , the 
equality holding only for pure initial states, p m = \ip)(ip\. 

Given the interferometer, i.e. given the operator J n to be measured, one can optimize 
the precision with which 9 is determined by choosing an initial state that maximizes the 
corresponding quantum Fisher information. 

In the case of distinguishable particles^ it has been shown that for any separable state 



7 Indeed, the inequality involving F[pi n ,J n ] given below (and similarly the set of inequalities discussed 
in |32j ) has been proved by using the decomposition of the collective operators J a , a — x,y,z, in terms 
of single-particle pseudo-spin: J a — Ja \ where refers to the z-th particle inside the trap. This is 
possible only if the particles are distinguishable. In the case of identical particles, single-particle operators 
like Ja" 1 are not addressable, otherwise one would be able to distinguish the particles by mean of them: only 
observables that are symmetric with respect to all permutations of the N bosons, like the J Q 's, are in fact 
permitted. 
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p sep the quantum Fisher information is bounded by N |25j : 

F[ Pscp ,J n ]<N . (26) 

This means that by feeding the interferometer with separable initial states, the best achiev- 
able precision in the determination of the phase shift 9 is bounded by the so-called shot- 
noise-limit (also called the standard-quantum-limit): 

A6 > — = . (27) 
~ y/N 



On the other hand, quite in general, one finds 

F[p,J n ]<N 2 , (28) 

so that an accuracy in phase estimation better than the shot-noise-limit is in principle 
allowed, eventually reaching the so-called Heisenberg limit, A6 > 1/N, obtained when in 
( 12 8p the equality holds. And indeed, many efforts have been devoted in order to find suitable 
input states p in and detection protocols that would allow such an ultimate sensitivity [7]- 
[55] . Notice that, because of f l2"B"|) . these states must be entangled. Actually, one can turn 
the argument around and use the quantum Fisher information for entanglement detection 
[25] ; indeed, if for a state p one finds that F[p, J n ] > N, than the state is surely entangled. 

When dealing with identical particles, as in the case of the condensed bosonic atoms con- 
fined in a double-well optical trap@ these conclusions need to be re-qualified. As discussed 
in the previous Sections, the notion of separability requires the choice of an algebraic bipar- 
tition; for the case at hand, the spatial bipartition (A±,A2) is a natural one: local operators 
are those that can be expressed as the product of polynomials in the creation and annihila- 
tion operators referring to the two wells. As shown earlier, with respect to this bipartition, 
a generic separable mixed state p is diagonal in the Fock basis ([7j), and can thus be written 
as in (Q. For such a state, the quantum Fisher information can be explicitly computed: 



TV N 

(29) 



F[ Pj J n ] = (n 2 x + n 2 y ) N + 2j2PkHN-k)-4Y, T* 1 (k + l)(N-k) 



In particular, in the case of a pure state, pk = \k,N — k)(k,N — k\, this expression results 
proportional to the variance of J n (cf. (|25l) ). 



F [p k , J n ] = (n 2 x + n 2 y )[N + 2k(N- k)] , (30) 

and can always be made greater than iV with a suitable choice of k. More specifically, 

i 2 



when n lays in the plane orthogonal to the z direction, so that n 2 + n 2 = 1, one finds 



8 When identical bosons are enough far away from each other so that their wavefunction do not spatially 
overlap, they effectively behave as distinguishable particles, and no symmetrization of the total wavefunction 
is needed 42 . However, this situation is hardly applicable to a gas of condensed ultracold atoms, where all 
particles share the same quantum state. 



10 



that F[pk, Jn\ > N for all k. Recalling flMj) . this implies that in this case the phase un- 
certainty is smaller than 1/yfN, thus beating the shot-noise-limit. Actually, when the 
two wells are filled by the same number of particles, so that the system is in the state 
Pn/2 = \N/2, N/2)(N/2, N/2\, one can even get close to the Heisenberg limit, since in this 
case: 

N 2 

F[p N/2 ,J n ]=— + N . (31) 

Therefore, unlike in the case of distinguishable particles, the quantum Fisher information 
can attain a value greater than N even with initial states that are separable with respect to 
the spatial bipartition. As a consequence, in general, the inequality (I2"fi| does not play any 
more the role of a separability condition when dealing with systems of identical particles. 

In spite of this, we have just seen that the accuracy with which the phase change can 
be determined in interferometers fed with such separable states can still beat the shot- 
noise-limit, provided that the rotation involved in the apparatus is not directed along the z 
axisjl As observed before, such rotations are realized by operators that are non-local with 
respect to the spatial bipartition0 Therefore, given the {A\,A2) bipartition, it is not the 
entanglement of the states fed into the interferometer that help overcoming the shot-noise- 
limit in the phase estimation accuracy; rather, it is the non-local character of the rotations 
operated by the apparatus on initially separable states that allows A8 to be larger than 
1/yN, with the possibility of closely approaching the Heisenberg 1/N limit. 

This result can be physically interpreted in another, equivalent way, which will shed further 
light on the notion of separability when dealing with identical particles. The idea is to change 
description (and thus bipartition) through a suitable Bogolubov transformation, following 
the discussion at the end of the previous Section. 

Take the unit vector n to lay in the plane orthogonal to the z axis, so that one can write 
n = (cosy?, simp, 0), p> G [0,27r]. Then, the generator J n in (l2"3"j) assumes the form: 

Jn = l(e- i va\a 2 + e i *a 1 a£) , (32) 

which clearly shows that its exponential is non-local in the {A\, A 2 ) bipartition. Nevertheless, 
it can become local in a different, suitably chosen bipartition. To this aim, let us introduce 
a new set of mode operators b\, bi, i — 1, 2 through the following Bogolubov transformation, 
that slightly generalizes the one in (Ql 



ai + e lip a 2 a x -e ^a 2 . . 

b ' = V2 • h= V2 ' <33) 



9 When n = (0,0,1), the quantum Fisher information (|29p vanishes. From the physical point of view, 
this results follows from the fact that the Fock states in (JT]) are eigenstates of the operator J z , so that 
separability is preserved by rotations generated by it: e ieJz \k,N — k) = e ie ( 2k ~ N )\k, N — k). In other words, 
in order to take advantage of the improvement in the accuracy of the phase determination, one has to use 
an experimental setup for which n ^ (0, 0, 1). 

10 In fact, only the exponential of J z happens to be a {Ay, ^2)-local operator, as shown in (TIT)) . 
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together with their hermitian conjugates. In this new representation, one has: 

Jn= l -{b\b l -b\b 2 ) , (34) 

so that the unitary operator that implements the rotation around n, 

e i0J n = e m\b x ii e -mlb 2 /2 ^ e ^ ^ e -tf&Jk/2 e ^ ^ ^5) 

is indeed local with respect to the new bipartition (81,82), where 8\ is the subalgebra of 
polynomials in b\, b\, while 82 is the one of polynomials in b\ , 62- I n this new language, 
the state |iV/2, N/2) representing the situation of equal filling of the two wells, is no longer 
separable with respect to this new bipartition; in fact, one finds (cf. (|2~Tj) ): 

W2.AV2) = E (f) (^f )(-!)"/- (4) N -*-'|0) , (36) 

while, as seen in the previous Section, any pure (B±, £>2)-separable state must be a Fock state 
of the form b\ b\ |0). 

Despite these changes, the value of the quantum Fisher information for the initial state 
\N/2,N/2) and the observable J n is unchanged and still given by fl3Tj) . since it does not 
depend on the representation used to compute it. This means that if one is able to build an 
experimental setup, together with a suitable measure procedure, which can be modelled in 
terms of the "energy" modes b\, bi instead of the original "spatial" modes a\, cii, then the 
accuracy A6 with which the phase 9 may be determined can still approach the Heisenberg 
limit. In such a case, the improvement in sensitivity with respect to the standard shot- 
noise-limit is due to the [B\, i32)-entanglement of the initial state \N/2, N/2) and not to the 
non-locality of the transformation that takes place inside the apparatus. 



5 Discussion 

The standard notion of separability for a many-body system made of N distinguishable 
particles is based on the natural tensor product structure of the Hilbert space in terms of 
the single-particle Hilbert spaces: a state (density matrix) of the iV-body system is separable 
if it can be written as the convex combinations of products of single-particle density matrices. 

In the case of a system of identical particles, due to the symmetrization (or antisym- 
metrization, for fermions) principle, this definition becomes meaningless. It can be replaced 
by a generalized one, that makes use of a "dual" language, focusing on the algebra A of 
operators of the system instead of the set of its quantum states. One fixes a partition of A 
in terms of a set of commuting subalgebras and defines as separable those states for which 
the associated expectation values of any factorized element of this partition can be written 
as a convex combination of products of expectation values. The notion of separability is thus 
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linked to a specific partition of A, so that a given many-body state can be separable with 
respect to one partition, but result entangled with respect to a different one. Nevertheless, 
this generalized definition of separability reduces to the familiar one expressed in terms of the 
single-particle tensor product structure in the case of a system of distinguishable particles. 

We have applied these considerations to the specific case of a system of N ultracold atoms 
trapped in an optical double-well potential, whose dynamics is very well captured by a two- 
mode Bose-Hubbard Hamiltonian. As we have seen, the second quantized language makes 
the application to this case of the new, generalized notion of separability very transparent 
and further allows the discussion of various related issues in quantum metrology. In fact, 
through state preparation and trapping potential control, this system has been shown to 
realize a highly sensitive Mach-Zhender interferometer, able to measure phase differences 
with a very high accuracy. 

Quite in general, the square error (A8) 2 in the determination of the phase difference 9 
accumulated inside the interferometer is bounded by the inverse of the quantum Fisher 
information F, whose value can not exceed N 2 . This gives the smallest possible error in 
the estimation of the phase, A9 > 1/N, the Heisenberg limit, which, for large N, is a huge 
improvement with respect to the standard shot-noise-limit, AO > l/y/N. 

In the case of a system of distinguishable particles, it has been proven that in order to 
beat the shot-noise-limit in the accuracy of the phase determination one needs to feed the 
interferometer with suitably iV-body entangled states. Indeed, one can show that for all 
separable states one has: F < N; as a consequence, the condition F > N signals the 
presence of entanglement and at the same time allows AO to be smaller than l/y/N^l 

When the interferometer is filled with identical particles, the condition F > N is no longer 
a univocal signal of state entanglement. Indeed, we have explicitly seen that in this case a 
quantum Fisher information larger than N may be obtained either via a non-local operation 
on separable states or via local operations on entangled states. Notice, however, that the 
notion of locality vs. non-locality and that of separability vs. entanglement need always to 
be referred to given algebraic bipartitions of the full algebra of observables. 

More specifically, in relation to the actual realization of a BEC-based interferometric ex- 
periment with a double- well optical trap, beam-splitting is usually implemented through 
lowering and raising of the inter-well potential barrier, while the subsequent phase estima- 
tion is obtained through number counting of particles inside the two wells. The algebraic 
bipartition that is relevant in this case is thus the spatial bipartition (Ai,A2), where local 

11 Many-body entanglement has also been related to spin squeezing [2Z]-[33]. For instance, in the case 
of distinguishable particles, a state for which the inequality N (AJ Z ) 2 > (J x ) 2 + (J y ) 2 is violated is surely 
entangled and spin squeezed, since the variance of J z is smaller than the standard quantum limit. Actually, 
a complete set of inequalities obeyed by all separable states have been discussed in [32] : violation of just one 
of them is enough to signal entanglement. Let us remark that in the case of a system of identical particles, 
some of these inequalities are violated even for separable states, thus loosing their role as entanglement 
witness. In this respect, when dealing with systems of identical particles, spin squeezing does not seem an 
unambiguous, useful resource for quantum metrology. A detailed discussion of these points will be reported 
elsewhere. 
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observables are those that can be expressed as the product of operators pertaining to the first 
and second well, respectively. A sub-shot-noise accuracy in the determination of the phase 
difference can then be obtained by acting with a non-local operation, i. e. a transformation 
generated by J x or J y , on mode-separable states, e.g. the balanced Fock state \N/2, N/2). 

Nevertheless, as shown at the end of Section 4, a different point of view can be equivalently 
adopted: it is based on an alternate measurement protocol, in which the energy bipartition 
(£>i, B2) becomes relevant. In this case, a local operation suffices to get a sub-shot-noise phase 
estimation accuracy, provided it acts on an entangled initial state. The practical realization of 
such a new type of BEC-based interferometer is surely an interesting experimental challenge. 
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